Control of recollision wave packets for molecular 
orbital tomography using short laser pulses 
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Abstract. The tomographic imaging of arbitrary molecular orbitals via high-order 
harmonic generation requires that electrons recoUide from one direction only. Within 
a semi-classical model, we show that extremely short phase-stabilized laser pulses offer 
control over the momentum distribution of the returning electrons. By adjusting the 
carrier-envelope phase, recollisions can be forced to occur from mainly one side, while 
retaining a broad energy spectrum. The signatures of the semi-classical distributions 
are observed in harmonic spectra obtained by numerical solution of the time-dependent 
Schrodinger equation. 
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A system driven by a strong laser field emits high-frequency radiation in a process 
called high-order harmonic generation (HHG). For atoms and molecules, HHG is well 
understood in terms of the three-step picture [1] as a sequence of field ionization, 
acceleration of the free electron in the laser field and recollision with the core. In recent 
years it has both become possible to use HHG to generate attosecond XUV pulses 
[21 [3l m [5l [6] , and to study molecular properties [3 [H [H [10]. Since the ground-breaking 
work of Itatani et al [11] it is known that HHG can be used to directly image molecular 
electronic orbitals in a scheme known as molecular orbital tomography. Itatani et al 
approximately reconstructed the highest occupied electronic orbital of N2 including its 
phase by comparing the high-harmonic spectra for different orientations of the molecule 
in the laser field with the spectrum of a reference system, namely atomic Ar, and 
combining this with knowledge about the electronic orbital of Ar. 

There are two main assumptions that enable us to perform the reconstruction of 
the molecular orbital; the first is that the returning electronic wave packets for the 
atomic and molecular case are very similar. This assumption stems from the fact that 
in the three-step picture the first step, tunnelling ionization, depends mainly on the 
instantaneous field strength and ionization energy. The two ionization energies are very 
similar. Also the second step, free electron acceleration, is identical for the two systems 
molecular N2 and atomic Ar if the binding potential is neglected. The second assumption 
is that the returning electrons in the third step of the three-step picture can be described 
by plane waves. This is supported by the fact that the classical trajectories leading to 
the higher harmonics involve spending roughly half a laser period in the continuum. 
This is enough time for the wave packet to spread out considerably to a diameter much 
larger than that of the core, leading to the effective plane-wave character. In reality, the 
incoming wave packet is distorted by the core potential, which will be a source of errors 
for the procedure. Due to the plane-wave approximation, the matrix element describing 
the recombination process leading to high-harmonic generation becomes proportional to 
a Fourier transform of the molecular orbital. The spectra for many different orientations 
of the molecule in the laser field and information about the continuum wave packet 
at the moment of recombination, obtained from the atomic reference system, can 
then be combined to reconstruct a projection of the molecular orbital onto the plane 
perpendicular to the the laser propagation axis. 

In general, recolliding electrons with both positive and negative momenta 
contribute to the generation of the same harmonic generation, each through their own 
recombination matrix element. For an arbitrary orbital, referring here to an orbital 
that does not have gerade or ungerade symmetry, there is no obvious way to deduct 
the matrix element needed for the reconstruction from knowledge about the harmonic 
spectrum and the continuum wave packet. The reconstruction of an arbitrary orbital 
therefore requires that the wave packets approach the core from one side only [121 fT3] . 
We propose achieving this behaviour of the continuum wave packet by using extremely 
short tailored laser pulses for high-harmonic generation. By stabilization and control 
of the carrier-envelope phase (CEP), i.e. the phase between the carrier wave and the 
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envelope of an extremely short pulse, the returning wave packet can be tailored. Because 
of the very-few-cycle character of the pulse, the CEP can be chosen such that classical 
trajectories in the unwanted direction are effectively suppressed. 

Return probability 

To characterize the returning wave packet we calculate semi-classically the probability 
that an electron returns with momentum k. We sample the ionization time with 3000 
points per laser cycle. At each ionization time ti, the electron tunnels with tunnelling 
probability (atomic units are used throughout) 



where E{ti) is the laser electric field at the moment of ionization and Jp is the ionization 
energy. After tunnelling, the electron follows a classical trajectory, starting with zero 
velocity at position zero. Every classical return receives a weight based on the tunnelling 
probability and a factor r"'^, where r is the time the electron spends in the continuum 
until the time of return. This factor reflects the effect of wave-packet spreading and 
comes from the Lewenstein model [Hj. We then add the contributions from the different 
trajectories to arrive at the probability that the electron returns to the core with a 
momentum k in the range k\ — A/c/2 < k < ki + Afc/2, where k\ is the momentum of the 
i'th momentum bin and A/c is the width of a bin. In contrast with the full Lewenstein 
model, interferences between the different paths leading to the same harmonic are not 
included in this model. This interference is irrelevant for the purpose of achieving one- 
sided electron returns. 

We consider two types of analytical representations for extremely short and linearly- 
polarized laser pulses. The first type are 3-cycle pulses with a sin^-envelope 



where Eq is the maximum field amplitude, T is the laser period and is the CEP. We 
choose a laser wavelength of 780 nm and an intensity of 5 • 10^^ W/cm^ for both pulses. 
The length of such a pulse in terms of the full width at half maximum of the intensity 
(FWHM) is 2.84 fs. The second type are pulses with a Gaussian envelope. These pulses 
are in principle infinitely long but for the semi-classical calculation we consider only the 
centre 4 cycles of the pulses 



These pulses have a FWHM of 2.50 fs. Both mentioned values of the FWHM are 
comparable to one laser cycle (2.6 fs). The choice of simple analytical expressions for 
the time-dependent field strength leads to small but unimportant violations of the rule 




(1) 




(2) 




(3) 



j E{t)dt = for the pulses with a Gaussian envelope. 
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Figure 1. (colour online). Returning electron momentum distributions for the 
3-cycle sin^-pulses. (a) CEP = rad, (b) (j) = 0.257r rad, (c) (p = O.SStt rad and 
(d) (f) = O.TStt rad. The distributions are shown on a linear scale (bars), and on a 
logarithmic scale (crosses). 

We set the ionization potential equal to /p = 30.2 eV, in accordance with the 
ionization potential of a frequently used 2D model system for [7j. In figures [T] and 
[2] the probability that an electron returns with momentum k is depicted for the two 
different pulse types and different CEPs. The bars reflect that probability on a linear 
scale, the crosses on a logarithmic scale. For each pulse the sum of the probabilities is 
normalized to 1. The insets show the laser pulse used to generate the data. Subjected to 
an infinitely long pulse of the mentioned intensity and wavelength, the fastest returning 
electrons would have an energy of 3.17 Up corresponding to ~ ±2.57 au, where 
Up = Eq/IAuj"^) is the ponderomotive potential [1]. The figures show that the CEP of an 
extremely short pulse has a very large influence on the characteristics of the returning 
wave packet, and that tuning of the CEP allows one to select the characteristics that 
one desires. This is also apparent in figure [3l which shows the fraction of the electrons 
returning with negative momentum as a function of the CEP for the two pulses. For 
= rad, the pulses have their main ionizing half-cycle in the positive direction, which 
means that the electrons leave towards the negative direction. As a consequence, most 
electron returns occur with positive momentum then. 

An optimal pulse for use with the molecular orbital tomography should meet two 
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Figure 2. (colour online). Returning electron momentum distributions for 
the Gaussian pulses. (a) CEP = rad, (b) (j) ~ Q.2tt rad, (c) (j) = O.Stt rad and 
(d) (f) — O.Ttt rad. The distributions are shown on a linear scale (bars), and on a 
logarithmic scale (crosses). 
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Figure 3. Fraction of electrons returning with negative momentum as a function of 
the CEP. Plus symbols correspond to the 3-cycle sin^-pulse, crosses to the Gaussian 
pulse. 



requirements: (i) all returning electrons should come from the same side, i.e. zero 
probability for electrons returning with the opposite sign of the momentum; (ii) the 
energy spectrum of the returning electrons should be broad, so that the Fourier 
transform of the orbital can be accurately extracted from the HHG spectra [12]. The 
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results show that a relatively large range of CEPs that are suitable for molecular orbital 
tomography can be found for both pulses. For both pulses the results are very similar. 
This indicates that many types of phase-stabilized extremely short laser pulses will be 
suitable for the procedure. These pulses should have an ionizing half-cycle where the 
amplitude of the electric field is significantly higher than everything earlier, and the 
following half-cycle in the opposite direction should be comparable to the ionizing half- 
cycle in strength, to force the electrons to return. Everything after the 'returning'-half- 
cycle should be significantly lower in amplitude, otherwise electrons will return from the 
other side with a significant probability. The reason that these effects can be achieved 
with the few-cycle pulses proposed in this paper, is that the ionization rate depends 
exponentially on the field amplitude. A nice demonstration that small field- amplitude 
changes can have large influence on harmonic spectra can be found in [15]. Because of 
this effect, we expect no problems with experimental pulses that have more smeared out 
tails, or with pre- and postpulses. For a different /p the qualitative behavior will be the 
same, as /p has no influence on the classical trajectories. Small differences in relative 
ionization rates might lead to a slightly different CEP being optimal for a different Jp 
though. 

Harmonic spectra 
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Figure 4. Harmonic spectra calculated using 3-cycle sin^-pulses with a CEP of 
(p = 0.257r rad (solid line) and </) = O.SStt rad (dashed line). The latter spectrum 
was multiplied with a factor 10 to separate the two curves graphically. 

A quantum-mechanical signature of the above-mentioned effects can be found in 
the high-harmonic spectra generated with these pulses. For this purpose we solve the 
time-dependent Schrodinger equation numerically using the split-operator method for a 
2D He+ model ion. We use the soft-core potential 
. 2 
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with the softness parameter adjusted such that Jp = 30.2 eV. In figure H] harmonic 
spectra associated with two values of the CEP for the 3-cycle sin^-pulse are shown. To 
ensure convergence, the wave function has been propagated for two more cycles after 
the end of the pulse. In both spectra, the low-frequency part is rather irregular due to 
the interference of several trajectories, while the high-frequency part shows the typical 
regular pattern of interference between short and long trajectories [16]. The latter is 
the pattern that we expect for uni-directional recollisions. 

For the case of = 0.257r rad, the transition to the regular interference pattern 
occurs at roughly the 28th harmonic. To relate the energy of an emitted photon Vt with 
the momentum of the returning electron /c, we use the energy-conserving dispersion 
relationship Vl = /c^/2 + /p. The 28th harmonic then corresponds to /c ~ ±1 au. As can 
be seen in figured! |A;| ~ 1 au is exactly the limit above which contributions come only 
from electrons with positive momentum. Also the cutoff at the 60th harmonic can be 
related to the classical picture, corresponding to the cutoff of /c ~ 2.2 au for electrons 
with positive momentum. The second cutoff comes from the less-likely events involving 
electrons with negative momentum, and extends clearly a couple of harmonics beyond 
the classically calculated and expected cutoff at the 76th harmonic, corresponding to 
k ~ —2.57 au. 

In the case of = O.SSvr rad, the classical picture suggests a single cutoff at 
k ~ 2.4 au, corresponding to the 68th harmonic. This is exactly what we see in figured] 
However, there is a small discrepancy between the TDSE and classical result concerning 
the location of the transition to a more regular spectrum. In the classical calculation, 
the fastest electrons with negative momentum have k ^ —1.3 au. This would correspond 
to a transition in the spectrum around the 33th harmonic. The TDSE result shows this 
transition around the 37th harmonic, which may be due to modifications of the electron 
trajectories by the Coulombic potential. 

Conclusion 

A necessary condition for accurate tomographic reconstruction of arbitrary molecular 
orbitals is that all recombining electrons return from the same direction. We have shown 
that this is achieved by using extremely short phase-stabilized laser pulses. Because of 
the exponential dependence of the tunnelling rate on the field amplitude, the choice of 
envelope shape is not very critical. 
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